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We study quantum degenerate ultra-cold Fermi gases in the presence of artificial spin-orbit cou- 
pling and crossed Zeeman fields. We emphasize the case where parity is violated in the excitation 
spectrum and compare it with the simpler situation where parity is preserved. We investigate in 
detail spectroscopic properties such as the excitation spectrum, the spectral function, momentum 
distribution and density of states for the cases where parity is preserved or violated. Similarly, we 
show that thermodynamic properties such as pressure, chemical potential, entropy, specific heat, 
isothermal compressibility and induced spin polarization become anisotropic as a function of Zee- 
man field components, when parity is violated. Lastly, we discuss the effects of interactions and 
present results for the pairing temperature as the precursor for the transition to a superfluid state. 
In particular, we find that the pairing temperature is dramatically reduced in the weak interaction 
regime as parity violation gets stronger, and that the momentum dependence of the order parameter 
for superfluidity violates parity when crossed Zeeman fields are present for finite spin-orbit coupling. 

PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 



I. INTRODUCTION 

The study of quantum degenerate fermions has been 
in the forefront of research in ultra-cold atoms and 
molecules in recent years, where particular attention was 
paid to the so-called evolution from BCS to BEC super- 
fluidity. Bringing ultra-cold fermions into quantum de- 
generacy and using Feshbach resonances to tune interac- 
tions between colliding fermions opened the door for the 
exploration of their superfluid phases, and their thermo- 
dynamic and correlation properties. This ability to tune 
interactions and explore the limits of weak, strong and 
unitary interactions had impact not only in cold-atom 
physics, but also in condensed matter physics (strongly 
correlated superconductors), nuclear physics (supercon- 
ductivity in quantum chromodynamics) and astrophysics 
(superfluidity in neutron stars). 

Many advances in cold atoms followed after the devel- 
opment of new tools for their toolbox. For instance the 
experimental study of the evolution from BCS to BEC 
superfluidity occurred after appropriate Feshbach reso- 
nances for 6 Li and 40 K were identified, and used to study 
the crossover problem for s-wave superfluids. In addition 
to s-wave Feshbach resonances, both 6 Li and 40 K also ex- 
hibit p-wave Feshbach resonances, which could produce 
p-wave superfluids in particular in BEC regime if p-wave 
Feshbach molecules were stable. Unfortunately, such p- 
wave molecules do not live sufficiently long for the cre- 
ation of p-wave superfluids (Tj-Q • 

In addition to the manipulation of interactions, it 
has been possible to extract experimentally detailed 
thermodynamic information of interacting ultra-cold 
fermions [5, 6] from local density images of trapped atoms 
with the help of the local density approximation 0, [1| . 

A further tool was developed recently through the pro- 
duction of tunable artificial spin-orbit fields [11] that were 
created in 87 Rb, a bosonic isotope of Rubidium, by using 
a set of Raman beams that allowed momentum transfer 



and mixing of two dressed spin states. In these experi- 
ments, the emergence of spin-orbit fields, controlled by 
the momentum transfer from the light fields to the atoms, 
was connected to the simultaneous existence of an artifi- 
cial Zeeman field controlled by the Raman coupling. In 
such experiments, the artificial spin-orbit and Zeeman 
fields were used to manipulate the effective interactions 
between bosons and thus produce new quantum phases 
of bosonic matter PJJ. In a subsequent experiment it 
was directly demonstrated that the effective interactions 
between bosons can acquire higher angular momentum 
components which are directly controlled by the artificial 
spin-orbit and Zeeman fields. The effective single-atom 
Hamiltonian created in these experiments is 

H sa (k) = e k l - h\\(j z - h x (k)a x - h y (k)a y , (1) 

when expressed in the dressed state spin basis |k, s). 
Here, ek = k 2 /(2m) represents the kinetic energy of an 
atom with momentum k and spin state s, hn = h z = Q/2 
represents the Zeeman field along the quantization axis 
z which is directly related to the Raman coupling SI, 
h x (k) = and /^(k) = h y + vk x , where h y = 5/2 repre- 
sents the detuning and vk x represents a mixture of equal 
Rashba [9[ and Dresselhaus [lOj terms, which we label as 
ERD spin-orbit coupling. 

The possibility of studying similar phenomena with 
ultra-cold fermions in the presence of spin-orbit cou- 
pling and Zeeman fields [HI, 03 lead to an explosion 
of theoretical research, which focused primarily on the 
case of Rashba-only (RO) spin-orbit coupling [l3l-[l7j 
which has been studied extensively in condensed matter 
physics in the context of non-centro-symmetric supercon- 
ductors [l8l-[20j . Some suggestions of possible realizations 
of artificial RO spin-orbit fields in ultra-cold fermions 
have appeared in the literature [21, 22]. However, in our 
group, we have focused mostly on the ERD case [23l — 
[26[, which is simpler to be realized experimentally as 
it was demonstrated for interacting bosons [U, [l2| in 
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the case of 87 Rb, and more recently for no n- interacting 
fermions such as 40 K in China [27] and 6 Li in the United 
States 

In our recent work [23l-[26|, we emphasized the impor- 
tance of interactions in producing novel phases for spin- 
orbit coupled ultra-cold fermions in the presence of Zee- 
man fields, and we analyzed in detail the ERD case with 
a longitudinal Zeeman field which is orthogonal to the 
spin-orbit field (23-l26j. In this paper, we generalize our 
previous analysis to include an additional Zeeman field 
h y representing the detuning J, which is perpendicular to 
longitudinal Zeeman field h\\ = h z representing the Ra- 
man coupling Q. In this case, the simultaneous presence 
of h y and the ERD spin-orbit field vk x along the y axis 
leads to the loss of parity for the eigenvalues of H sa (k) de- 
fined in Eq. (pQ) above. This violation of parity in the en- 
ergy spectrum leads to many detectable parity- violating 
properties in spectroscopic quantities such as the spectral 
function and momentum distribution, which are explic- 
itly momentum-dependent, or in momentum-integrated 
quantities like the density of states when viewed as a 
function of Zeeman fields h y and h y . Furthermore, we 
demonstrate that this parity violation can also be seen 
and quantified in thermodynamic quantities such as the 
pressure, entropy, specific heat, chemical potential, com- 
pressibility and spin-polarization, which are discussed in 
detail for non-interacting ultra-cold fermions in the pres- 
ence of ERD spin-orbit and crossed Zeeman fields h y and 
h z . Furthermore, we also discuss the effects of parity vi- 
olation when attractive s- wave interactions are present, 
and we pay particular attention to the pairing tempera- 
ture and to the absence of parity in the order parameter 
tensor describing the supernuid order. 

The remainder of the paper is as follows. In section HI} 
we describe the magnetic hamiltonian in its general form, 
and particularize it to any linear combination of Rashba 
and Dresselhaus terms, while in section HIJ we discuss the 
independent-atom (non-interacting) Hamiltonian includ- 
ing both the kinetic energy and the magnetic parts. The 
eigenvalues and eigenvectors of the single-atom Hamilto- 
nian are discussed in section IIV1 where a generalized he- 
licity basis is introduced and Fermi surfaces for various 
values of crossed Zeeman fields are presented. Particu- 
lar attention is paid to parity violations for the eigen- 
values and eigenvectors. In section [Vj we show the mo- 
mentum dependence for the spectral function at fixed 
frequency revealing the cases of weak and strong parity 
violations, while in section PVTI we describe the the spin- 
resolved momentum distributions which also reveal par- 
ity violations for various values of crossed Zeeman fields. 
In section IVII1 we analyze the spin-dependent density of 
states for a few values of the crossed Zeeman fields and 
fixed ERD spin-orbit coupling, while in section IVlIIl we 
define several thermodynamics properties which analyzed 
in the following sections, such as pressure and entropy in 
section IIX| chemical potential in section [Xj isothermal 
compressibility in section IXII and spin-polarization IXHI 
Furthermore, we investigate the effects of parity viola- 



tion on the effective interactions in section IXIIII and on 
the pairing temperature in section IXV[ Lastly, we state 
our conclusions in section IXVII emphasizing that several 
measurements can be made to detect parity violation in 
both non-interacting and interacting Fermi gases in the 
simultaneous presence of an ERD spin-orbit field, and 
crossed Zeeman fields h y , h z , where h y is related to the 
detuning 5 and h z to Raman coupling Q. 

II. MAGNETIC HAMILTONIAN 

To describe quantum degenerate Fermi systems in the 
presence of artificial spin-orbit and crossed artificial Zee- 
man fields, we discuss first the type of magnetic Hamilto- 
nian to be used. Generally speaking the coupling between 
magnetic fields and two spin-states is given by H mag (k) = 
— h(k) • <r, where h(k) is the effective magnetic field (in- 
cluding Zeeman and spin-orbit) and a = (a x ,a y ,a z ), 
with <Ji being the Pauli matrices. This leads to the mag- 
netic Hamiltonian 

rr (Xc \-( M k ) h x (k) - ih y (k)\ { , 

^ ^ " \h X (k) + ihy(k) ~h Z (k) J ' W 

Although there are many possible forms of fictitious 
magnetic fields that can be created in the laboratory. We 
discuss here a few simpler magnetic field configurations. 
In principle, the external artificial Zeeman field Icize can 
have three components (h x ,h y ,h z ), in practice, we set 
h x = 0, and (h y , h z ) ^ because in current experimental 
setups h y = 5/2, where 5 is the laser detuning, and h z = 
Q/2, where £1 is the strength of Raman coupling field. 

In addition, there are many possible types of spin- 
orbit contributions where hso can have three compo- 
nents (hso,x(k),hso,y(k),hso,z(k)). However, we con- 
sider particular forms of spin-orbit fields, which can be 
more easily realized in practice. The first type has the 
Dresselhaus form ho(k) = vd (k y , k x , 0) , where vd 
measures the strength of the Dresselhaus field in units of 
velocity. The corresponding Hamiltonian for such field is 

*^ = - v °{k y lik x ky 7 kx )- ^ 

The second type has the Rashba form h#(k) = 
VR^—ky, k x , 0) , where vr measures the strength of the 
Rashba field in units of velocity. The corresponding 
Rashba Hamiltonian is 

H *( k )=^GA* kv+ o lkx )- (4) 

Either the Dresselhaus or the Rashba forms require spin- 
orbit fields along the x and y directions, which to be 
produced experimentally demand two orthogonal Raman 
setups, such that momentum transfer occurs in two per- 
pendicular directions. A linear combination of two forms 
leads to the field h RD (k) = (v RD _k y , v RD+ k x , 0) , 
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where the velocities vrd± = (vd ±vr). The correspond- 
ing Hamiltonian for such linear combination is 

w v RD _k y +iv RD+ k x \ 

UrdW- [ VRD _ ky - ivRD+kx )' 

(5) 

The simplest type of spin-orbit field that can been cre- 
ated in the laboratory has the equal- Rashba-Dresselhaus 
(ERD) form h £ ; jRjD (k) = v(0,k x ,0), where v RD _ = 0, 
and vrd + = v, or equivalently vd = v R = v/2. The 
corresponding Hamiltonian for the ERD spin-orbit field 
has the simple form 

Taking into account an arbitrary superposition of 
Rashba and Dresselhaus spin-orbit coupling and a gen- 
eral uniform Zeeman field Icize = (h x ,h y ,h z ), we can 
write the Zeeman-spin-orbit Hamiltonian as 

where the parallel component of the total field is hu = h z 
and the transverse component is 

ft±(k) = [h x + h RDyX (k) - i(h y + h RDi y(k)] , 

where h RDjX (k) = v RD _k y and h RDjy (k) = v RD+ k x . 

Having presented the magnetic Hamiltonian, we dis- 
cuss next the Hamiltonian including the kinetic energy 
of the atoms. 



III. HAMILTONIAN 

The Hamiltonian for non-interacting ultra-cold 
fermions with identical masses m in the presence of 
spin-orbit and crossed Zeeman fields can be written in 
second quantization as 

W = ^*t( k ) H (k)*(k), (8) 

k 

where the spinor ^(k) = ^cj^cj^ describes the cre- 
ation of fermion states with momentum k and spin t or 
Such Hamiltonian describes two hyperfine states of 
Fermi atoms such as 6 Li or 40 K, and the corresponding 
Hamiltonian matrix is 

where K a (k) = k 2 /(2m) — represents the kinetic en- 
ergy of a fermion with mass m, momentum k and spin a 
with respect to the chemical potential fi a . 

We define the variables K± = (^K^ ± K^j /2 and the 

chemical potentials fi± = (/x^ ± m) /2, and notice that 



K + = |k| 2 / (2m) — [i+ = plays the role of the average 
kinetic energy, while K- = — + h z ) = — fty, plays the 
role of the parallel Zeeman field hu including the external 
field h z and the internal field \i- due to a possible initial 
population imbalance. In the limit that there is zero 
initial population imbalance, the chemical potentials can 
be set to /i + fi and /i_ — >• 0, while the kinetic energies 
reduce to K+(k) — >> £k and K-(k) — >> —h\\ = —h z , where 
£ k = |k| 2 /(2m) - [i. 

The non-interacting Hamiltonian has the simpler form 
Ho(k) = ^kl — h\\a z — h x (k)a x — h y (k)a y which can be 
re-expressed as 

where h±(k) = h x (k) — ih y (k). In the remainder of 
the manuscript, we will discuss this explicit form of the 
Hamiltonian and some interesting consequences. 

We define the total number of fermions as TV = 
iV-f + iVj,, and choose our energy, velocity and momentum 
scales through the Fermi momentum kp defined from the 
total density of fermions n = + = /4/(3tt 2 ), where 
n = N/V and n s = N s /V with s = (t,|). This choice 
leads to the Fermi energy ep = k F /2m and to the Fermi 
velocity vp = kp/m, which are the energy and velocity 
scales used throughout the manuscript. 

IV. EXCITATION SPECTRUM 

Now, let us introduce the unitary matrix Uk that di- 
agonalizes the Hamiltonian Ho(k), such that 

E(k) = u[H (k)U k (11) 

is a diagonal matrix containing the eigenvalues of Ho(k). 
The corresponding eigenvectors are the spinors 3>(k) = 
Uj c ^ r (k), where the unitary matrix Uk has a momentum- 
dependent SU(2) form and can be written as 

°- =(-<£)• <12 » 

where the normalization condition \u\r\ 2 + |^k| 2 = 1 is 
imposed to satisfy the unitarity condition U^Uk = 1 
leading to the following expressions 

where is taken to be real without loss of generality 
and 

is a complex function where the phase cp^ is defined via 
h ± (k) = |/i ± (k)|e^ k , leading to <p k = Arg [h ± (k)} . 



The eigenvalues of H (k) emerge as 

E « = C*o k) J k) ). (») 

where £ft(k) = ^ — |h e ff(k)| is the eigenvalue where the 
momentum dependent effective field h e ff(k) is aligned 
with the spin ff, and ^(k) = £ k + |h e ff(k)| is the 
eigenvalue where the momentum dependent effective field 
h e ff(k) is aligned with the spin 4- Here, |h e ff(k)| = 

^h 2 + |/i^(k)| 2 is the magnitude of the effective field. 

The respective eigenvectors are 

<^(k) = ^ k c kt - v k c k ; (16) 
corresponding to the state |k, ff), and 

$ 4 (k) =<c kt + w k c k; (17) 
corresponding to the state |k, J|). 
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FIG. 1. (color online) Helicity energy dispersions e^(k)/eF 
(blue line) and e^(k)/eF (magenta line) versus momentum 
k x /kF with k y — k z — for ERD spin-orbit coupling v/vf — 
0.71. For reference, the black dashed lines show the helicity 
bands for v/vf = 0.71 with h z /eF = h y /cF = 0. The Zeeman 
fields are (a) h y /cF = and h z /eF = 0.1, (b) h y /cF — and 
h z /eF — 1.2, (c) hy/eF = 0.1 and h z /eF = 0, (d) h y /tF — 0.5 
and h z /eF = 0, (e) h y /cF = h z /eF = 0.1, (f) h y /cF = 0.5 
and h z /eF — 1.2. 

For initially population balanced systems in the ERD 
case, with the Zeeman field having only and h z com- 
ponents, the magnitude of the effective magnetic field 

\h G fi (k)| = J h 2 + (/fy + v/c^) 2 does not have well defined 



parity, which implies that the same is true for the eigen- 
values £a(k). If either h y = (zero detuning) or v = 
(zero spin-orbit coupling) , parity is restored for the eigen- 
values. This can be seen in Fig. [TJ where the energy dis- 
persions e^(k) = €k - |ft e ff (k)| and e^(k) = e k + |ft e ff(k)| 
are shown for several cases. In all panels of Fig. [T] the 
dashed lines represent the case where the ERD spin or- 
bit is finite (v/vp = 0.71), but the Zeeman fields are 
zero (h z = h y = 0). This case corresponds to shifted 
parabolic bands e^(k) = k 2 /(2m) — \vk x \ + k 2 /(2m) 
and ejL(k) = k 2 x /(2m) + |^| + k 2 p /(2m), where ^ = 
k 2 + k 2 . The lower helicity band e^(k) has two min- 
ima, which occur at (fc^ = m\v\,k y = 0,k z = 0) and 
(fcar = — m\v |, fc^ = 0, /c 2 =0), respectively, while the up- 
per helicity band e^(k) has only one minimum occurring 
at (k x = 0, k y = 0, k z = 0). This can be seen by complet- 
ing the squares and rewriting the dispersions of the helic- 
ity bands as e^(k) = (\k x \ — m\v\) 2 / (2m) — e v + k 2 /(2m) 
and e^(k) = (|^| + m\v\) 2 /(2m) - e v + k 2 /(2m), where 
e v = mv 2 /2 is the characteristic kinetic energy associ- 
ated with the spin-orbit coupling strength v, which has 
units of velocity. 

In Fig. UK, we show the case of v/vp — 0.71, for h y = 
and h z /cF = 0.1, and in Fig. [Ud, we show the case of 
v/vf = 0.71, for h y = and h z /cF = 1.2. In both of 
these cases the energy dispersions are parity preserving, 
and the main fundamental difference between Fig. [T^i and 
Fig- lb is that the lower helicity band e^(k) has two min- 
ima at finite k x = \J ^e 2 — h 2 as seen in Fig. [TJi, 
but a single minimum at k x = as seen in Fig.flJ}. While 
the upper helicity band e^(k) has only a single minimum 
at k x = in both Fig. QJi and Fig. [TJd. The two min- 
ima in the lower helicity band for Fig. [TJi occur only for 
low Zeeman fields \h z \ < 2e v . In Figs. QJ and [T]i, we 
show the helicity bands for h z = 0, but h y /cF = 0.1 and 
hy/eF = 0.5, respectively. In these cases, the energy dis- 
persions are (k) = k 2 / (2m) — | h y + vk x | + k 2 / (2m) and 
ej| (k) = k 2 1 (2m) + \h y + vk x \ J rk 2 / (2m) and do not even 
parity as it is standard, since e a (k) ^ e a (— k). These dis- 
persions are compared to the dispersions for the case of 
hy/ep = h z /ep = 0, but finite v/vp, which are shown as 
dashed black lines in Fig. [U The last examples described 
in Fig. []J and Fig. [TJ correspond to cases where both h y 
and h z are non-zero, having values h y /ep = h z /ep = 0.1 
for (e) and h y /ep = 0.5, h z /ep = 1.2 for (f). Notice 
the presence of two minima for the lower helicity band in 
(e), and the existence of only one minimum for the lower 
helicity band in (f), but in both cases parity (inversion 
symmetry) is violated. 

In Fig. [2j we show cross sections of the Fermi surfaces 
(FS) in the k x -k y plane at k z = for ERD spin-orbit 
coupling and a few values of crossed Zeeman fields. The 
Fermi surfaces have rotational symmetry about the k x 
axis, that is, in the k y -k z plane, and their full three- 
dimensional structure can be visualized using this prop- 
erty. We show in Figs. [2ji-c the case for Zeeman fields 
hy/ep = and h z /ep = 0.1, where the Fermi surfaces 
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FIG. 2. (color online) Cross sections of Fermi surfaces 
in the k x -k y plane at k z — for ERD spin-orbit coupling 
v/vf = 0.71. In (a) through (c), the Zeeman fields are 
hy/eF = and h z /eF = 0.1, where parity symmetry is pre- 
served, and the different values of the chemical potential and 
induced polarization are (a) /jl/ef = 0.495, Pi n d = 0.121, (b) 
/jl/ef = 0, Pind = 0.147, (c) n/e F = -0.370, P ind = 0.106. In 
(d) through (f) the Zeeman fields are h y /cF = h z /eF = 0.1, 
where parity symmetry is not preserved, and the different val- 
ues of the chemical potential and induced polarization are (d) 
l^/ep = 0.492, P ind = 0.121, (e) /i/e f = 0, P ind = 0.143, (f) 
/jl/ef = —0.370, Pind = 0.100. The blue and magenta lines in- 
dicate the energy contours of the lower e^(k) and upper e^(k) 
helicity bands, respectively. 



exhibit parity (or inversion) symmetry. These parame- 
ters correspond to the helicity bands shown in Fig. QJi, 
where the lower helicity band has two minima. While we 
show in Figs. EJi-f the case for Zeeman fields h y /ep = 
h z /ep = 0.1, where the corresponding Fermi surfaces do 
not have well defined parity or inversion symmetry. The 
specific values of the chemical potential /i and induced 
polarization P[ n ^ = (N^ — N_i)/(N^ + N±) are indicated 
in the captions. 

In Fig. [2j notice also that as the chemical potential 
ji is changed below the bottom of the helicity band 
e^(k), the central pocket of the Fermi surface disappears 
(see magenta surface near zero momentum in Fig. [2]). 
Further lowering of the chemical potential leads to the 
crossing of a local maximum of the helicity band e^(k), 
where the residual Fermi surface break into two pock- 
ets (see blue surfaces in Fig. [2]). This is reminiscent 
of the Lifshitz transition [29j in non- interacting met- 
als, also called metal-to-metal or conductor-to-conductor 
transition, where under pressure or another external pa- 
rameter the Fermi surface of the system changes topol- 
ogy producing a major rearrangement of momentum 
states which lead to a drastic change in the density of 
states of the system. The thermodynamic potential fl 
in the vicinity of the usual Lifshitz transition behaves as 



Q = Q reg -f a\ji — /i c \ 5 ^ 2 , where ^ reg is the regular (an- 
alytic) part, and a is the prefactor of the non- analytic 
component. The isothermal compressibility is related 
to the second-derivative of the thermodynamic potential 
with respect to the chemical potential and behaves as 
kt = ^T,reg + — H^ 1 ^ 2 •> where ftr,reg is the regular 
part and f3 is the coefficient of the non-analytic com- 
ponent. According to Ehrenfest's classification of phase 
transitions, the non-analyticity manifests itself only in 
the third derivative and is a third-order phase transi- 
tion. However, the Lifshitz transition is more commonly 
called the 2-1/2 order transition in allusion to the spe- 
cific 5/2 power- law non-analyticity of Q in three dimen- 
sions. This topological transition is not characterized 
under Landau's symmetry-based classification, since no 
symmetry is broken in the Lifshitz case. This trivial Lif- 
shitz transition can be seen in Fig. [2] both for the parity- 
preserving and parity- violating examples. 



In Fig. [3j we also show cross sections of the Fermi sur- 
faces (FS) in the k x -k y plane at k z = for ERD spin-orbit 
coupling v/vp = 0.71 and a few values of crossed Zeeman 
fields. The Fermi surfaces have also rotational symme- 
try about the k x axis, that is, in the k y -k z plane. We 
show in Figs. [3^i-b the case for Zeeman fields h y /ep = 
and h z /ep = 1.2, where the Fermi surfaces exhibit parity 
(or inversion) symmetry. These parameters correspond 
to the helicity bands shown in Fig. [TJd, where the lower 
helicity band has only one minimum. While we show in 
Figs. [3]>d the cases of the Zeeman fields h y /ep = 0.5 
and h z /cp = 1.2 where the corresponding Fermi sur- 
faces do not have well defined parity or inversion sym- 
metry. The values of the chemical potential \i and in- 
duced polarization Pi n d are /i/ep = 1.5, Pi n d = 0.692 for 
Fig. Hi; n/e F = 0, P ind = 0.801 for Fig. Eb; /i/e F = 1.5, 
P ind = 0.673 for Fig. [3t; and /i/e F = 0, P ind = 0.739 
for Fig. [3jd. Notice that there is a fundamental difference 
between the Fermi surfaces in Figs. [2] and [3] in connec- 
tion with their topology. The lower helicity band e^(k) 
can have two simply connected FS for the parameters of 
Fig. [21 but only one simply connected FS for the param- 
eters of Fig. [3l This means that there is only one triv- 
ial Lifshitz transition in the case of Fig. [3j while there 
are two trivial Lifshitz transitions in the case of Fig. [2j 
We call this transition for non-interacting systems trivial 
to contrast it with a more exotic, but related topolog- 
ical transition that can occur in p-wave [HO, HH or d- 
wave [32|, [33[ superfluids, where interactions play a fun- 
damental role. 



A further characterization of the parity violation 
present in fermion systems with spin-orbit coupling and 
crossed Zeeman fields can be made by analyzing addi- 
tional spectroscopic quantities such as the spectral func- 
tion to be discussed next. 
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FIG. 3. (color online) Cross sections of Fermi surfaces in the 
kx-ky plane at k z = for ERD spin-orbit coupling v/vf = 
0.71. In (a) and (b), the Zeeman fields are h y /cF = and 
h z /eF = 1.2, with parity symmetry being preserved. The 
values of the chemical potential and induced polarization are 
(a) fi/e F = 1.5, Find = 0.692; (b) fi/e F = 0, P ind = 0.801. In 
(c) and (d), the Zeeman fields are h y /cF = 0.5, h z /cF = 1.2 
and parity symmetry is not preserved. The values of the 
chemical potential and induced polarization are (c) /jl/ef = 
1.5, P ind = 0.673; (d) fi/e F = 0, P ind = 0.739. The blue and 
magenta lines indicate the energy contours of the lower (k) 
and upper e^(k) helicity bands, respectively. 



V. SPECTRAL FUNCTION 

A very useful tool to probe momentum-resolved prop- 
erties is the use of radio-frequency (RF) spectroscopy 
that can extract the spectral function [34] yielding simi- 
lar measurements to those encountered in photoemission 
spectroscopy of condensed matter physics. The resolvent 
(or Green) operator matrix is defined as 



G(k,z) 



1 



zl - H(k) 



(18) 



in momentum- frequency space. In the present case, the 
diagonal components of G(k, z) are 



G t (M) 



z - £ t (k) 
for the spin j" component and 



^(k) 



Gi(k,z) 



Ct(k) 



(19) 



(20) 



for the spin 4- component. The corresponding spectral 
function is A s (k, uo) = — 7r _1 Im [G s (k, a; + iS)] , which in 



terms of the coherence factor and v k becomes 

4 t (k, u) = v£6(w - ^(k)) + \v k \ 2 5(u; - ^(k)) (21) 

for the up-spin (t) component, and 

^(k, u) = \v k \ 2 5(u; - ^(k)) + v£6(w - ^(k)) (22) 

for the down-spin (|) component. 

The spectral functions A s (k, gj) at frequency uo = n 
and temperature T/e^ = 0.05 are shown in Fig. |4] for 
some values of the crossed Zeeman fields and spin-orbit 
coupling. A small energy broadening 5/ep = 0.01 is in- 
cluded and a logarithmic scale is used to help visualiza- 
tion. In the relevant panels of Fig. HJ the blue dashed 
lines represent A s (k x = 0,k y ,k z = 0,cj = /i), and the 
red solid lines represent A s (k Xl k y = 0, k z = 0, u = /i). 
Additionally, the two left-most columns correspond to 
the Af(k) component, and the two right-most columns 
describe the Aj,(k) component. In Figs. H^a)-(d) the 
Zeeman fields are h y /cF = h z /eF = 0, the ERD cou- 
pling is v/vf = 0, the chemical potential is fi/ep = 1.0, 
and the induced polarization = 0. In Figs. H{e)- 
(h) the Zeeman fields are hy/cp = 0, h z /ep = 0.1, the 
ERD coupling is v/vp = 0.71, and the chemical poten- 
tial is fjb/ep = 0.495, and the induced polarization is 
^Rnd = 0.121. In these two cases parity is preserved, since 
the eigenvalues e a (k) and coherence factors and are 
invariant under momentum inversion for h y /ep = 0. In 
Figs.|3(i)-(1) the Zeeman fields are h y /ep = h z /ep = 0.1, 
the ERD coupling is v/vp = 0.71, the chemical poten- 
tial is ji/ep = 0.492, and the induced polarization is 
^Pind = 0.120. Parity is not preserved in the last case 
since (h y /ep ^ 0). This is reflected in the absence of 
inversion symmetry, which is noticeable but weak. 

In Fig. El we show the spectral functions A s (k, uo) at 
frequency uj = /i and temperature T/ep = 0.05 for ERD 
spin-orbit coupling v/vp = 0.71. For varying chemical 
potentials, we choose the particular values of the crossed 
Zeeman fields to be h y /ep = h z /ep = 0.1 in order to 
emphasize the absence of inversion symmetry (parity) 
when h y ^ 0. As before, a small energy broadening 
5 / 6p = 0.01 is included and a logarithmic scale is used 
to help visualization. The same color convention is used 
in the relevant panels of Fig. [5j where the blue dashed 
lines represent A s (k x = 0, k y , k z = 0, uj = /x), and the red 
solid lines represent A s (k Xl k y = 0, k z = 0, uj = /x). In 
addition, the two left-most columns describe the Af(k) 
component, and the two right-most columns correspond 
to the Aj^k) component. In Figs. [5fa) -(d) the chemical 
potential is \ijep = 0.492, the induced polarization is 
-Rnd = 0.120, and parity is weakly violated, but notice- 
able. In Figs.[5fe)-(h) the chemical potential is fi/ep = 0, 
the induced polarization is Pi n d = 0.143, and parity is 
more strongly violated. In Figs. Efi)-(l) the chemical 
potential is ji/ep = —0.37, the induced polarization is 
^ind = 0.105, and parity is strongly violated. 

In the last two cases parity is violated more strongly, 
since the eigenvalues e a (k) and coherence factors 
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FIG. 4. (color online) Finite temperature (T/cf = 0.05) 
dimensionless spectral functions 27V€FA s (k, uj) at uj = ji for 
parameters (a)-(d) v/vf = 0, h y /cF = h z /eF = and /h/cf = 
1.0; (e)-(h) v/vf = 0.71, h y /cF = 0, h z /eF = 0.1, and /x/e^ = 
0.495; (i)-(l) v/v F = 0.71, h y /e F = h z /e F = 0.1 and fi/e F = 
0.492. The blue dashed lines are cuts along the direction 
(k x = 0, k y , k z = 0) corresponding to 2neFA s (k x = 0,k y ,k z = 
0,cj = /i), and the red solid lines are cuts along the direction 
(k x , k y = 0, k z = 0) corresponding to 27Y€FA s (k x , k y = 0, k z = 
0,cj = /i). The two left-most columns describe the Af(k) 
component, and the two right-most columns correspond to 
the A^(k) component. 

and are not invariant under momentum inversion for 
hy/cp ^ and are more sensitive to this violation for 
chemical potentials closer to the bottom of the helicity 
bands. 



VI. MOMENTUM DISTRIBUTION 

To understand the momentum distribution for quan- 
tum degenerate Fermi gases in the presence of spin-orbit 
and crossed Zeeman fields, we look first at the expec- 
tation value (n a (k)) of the number operators n a (k) = 
(k)<I> a (k), which describes the momentum distribu- 
tion n a (k) in the helicity basis. In this case, the momen- 
tum distribution is n a (k) = np(£ a (k)), where uf(x) is 
the Fermi function. 

The momentum distributions for the original spin 
states s = (t?i) are defined as n s (k) = (cj^Cks). With 
the help of the unitary matrix U(k) defined in Eq. ([T2]h 
which relates the creation and annihilation operators in 
the helicity and the standard spin basis the momentum 
distributions become 

n t (k) = t£„ F (^(k)) + M 2 n F (a(k)) (23) 

for the up-spin (t) component, and 

% (k) = | Uk | V(^(k)) + u£M&(k)) (24) 




ki/k F kx/k F kx/k F 

FIG. 5. (color online) Finite temperature (T/cf = 0.05) 
dimensionless spectral functions 27V€FA s (k, uj) at uj = ji for 
ERD spin-orbit coupling v/vf — 0.71 and Zeeman fields 
hy/eF = h z /e F = 0.1 In panels (a)-(d) \i/e F = 0.492; (e)- 
(h) n/e F = 0; and (i)-(l) fi/e F = -0.37. The blue dashed 
lines are cuts along the direction (k x — 0,k y ,k z — 0) corre- 
sponding to 27V€FA s (k x = 0,k y ,k z = 0,cj = /i), and the red 
solid lines are cuts along the direction {k Xl k y — 0, k z — 0) 
corresponding to 2neFA s (k X) k y = 0, k z = 0, cj = ji). The two 
left-most columns describe the A^(k) component, and the two 
right-most columns correspond to the A|,(k) component. 

for the down-spin (\) component. 

Such expressions can be also obtained from the general 
relation 

/oo 
dujnF(oj)A s (\<L, u), (25) 
-oo 

between the momentum distribution and the spectral 
function for fermions. 

It is also convenient to obtain the momentum distri- 
bution sum n+(k) = n^(k) + n^(k) and the momentum 
distribution difference n_(k) = n^(k) — n^(k). The first 
distribution can be written in terms of the Fermi func- 
tions only 

n+(k) = n F (^(k)) + n F (^(k)), (26) 

while the second distribution can be written in terms of 
the Fermi functions and the components of the effective 
Zeeman field 

n-(k) = K(^(k) - n F (&(k)] . (27) 

In Fig. [6j we show momentum distributions n s (k) at 
T/cf = 0.05 and in the regime where the Fermi system is 
largely degenerate, containing wide regions in momentum 
space where n s (k) « 1. The blue dashed lines represent 
cuts of n s (k) along the k = (k x = 0,k y ,k z = 0) direc- 
tion, while the red solid lines represents cuts of n s (k) 
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along k = (k xi k y = : k z = 0). The left-most columns 
correspond to spin ^, while the right-most columns rep- 
resent spin |. For reference, we show in Figs.[6ja)-(d) the 
case with zero spin-orbit coupling and without Zeeman 
fields, corresponding to parameters v/vf = 0, h y /eF = 
h z /cF = chemical potential fi = 1.0. In this case, the 
momentum distributions for the two spin components 
are identical n^(k) = n^(k), meaning that the popula- 
tions are balanced with n_(k) = 0, and induced polariza- 
tion is Pi n d = 0. In Figs. [6je)-(h), we show momentum 
distributions for parameters v/vp = 0.71, h y /ep = 0, 
h z /ep = 0.1 and chemical potential \i = 0.495. Here, the 
momentum distributions acquire double plateaux struc- 
tures along the direction k = (k x , k y = 0, k z = 0) due to 
the momentum shifts of the helicity bands as shown in 
Fig. [Ha). Additionally, the momentum distributions for 
different spin-components are no longer identical, such 
that n^(k) ^ ^H(k), or n_(k) ^ 0, and the induced po- 
larization is non-zero taking the value Pi n d = 0.121. In 
these two cases, parity is not violated and the momentum 
distributions are even functions of momentum under in- 
version symmetry. However, in Figs. EJi)-(l), we show 
momentum distributions for parameters v/vp = 0.71, 
hy/eF = h z /cF = 0.1, and chemical potential \i = 0.492, 
in which case the double plateaux structures are still pre- 
served, population imbalance is present with n_(k) ^ 
and induced polarization P in & = 0.120. Most importantly 
parity is weakly violated since n s (— k) ^ n s (k). 

In Fig. [7J we show momentum distributions n s (k) at 
T/cf = 0.05 for parameters v/vp = 0.71, h y /ep = 
h z /ep — 0.1, and varying chemical potentials p/ep. We 
emphasize the regimes where parity is more strongly vio- 
lated leading to momentum distributions without inver- 
sion symmetry: n s (— k) ^ n s (k). The blue dashed lines 
represent cuts of n s (k) along the k = (k x = 0, k yi k z =0) 
direction, while the red solid lines represents cuts of 
n s (k) along k = (k Xl k y = 0, k z = 0). For reference, 
we show in Figs. [T^a)-(d) the case with ji/ep = 0.492 
and Pi n d = 0.120, where the Fermi system is still largely 
degenerate, containing wide regions in momentum space 
with n s (k) w 1, and at the same time parity is vio- 
lated only weakly. In Figs. [T^e)-(h), we show momen- 
tum distributions for fi/cp = and Pi n d = 0.143, while 
in Figs. [Tfi)-(l), we show momentum distributions for 
li/ep = —0.37 and Pi n d = 0.105. In both cases, the 
Fermi system is no longer degenerate, containing wide 
regions in momentum space where n s (k) < 1. In the 
last two cases the momentum distributions remain sym- 
metric upon reflection along the k y or k z directions, but 
parity is more strongly violated leading to a highly asym- 
metric momentum distributions along the k x direction. 

Having discussed the momentum distribution at low 
temperatures, we analyze next the density of states of 
quantum degenerate fermions in the presence of spin- 
orbit coupling and crossed Zeeman fields. 




ki/k F kx/k F kx/k F ki/k F 



FIG. 6. (color online) Finite temperature (T/cf — 0.05) 
momentum distributions n^(k) (two left-most columns) and 
ni(k) (two right-most columns). The parameters for (a)-(d) 
are ERD spin-orbit coupling v/vf = 0, Zeeman fields h y /cF = 
h z /eF — 0, chemical potential /j/cf = 1.0, and induced po- 
larization Pind = 0. Similarly, the parameters for (e)-(h) are 
v/vf — 0.71, hy/e = 0, h z /eF = 0.1, h/ef — 0.495, and 
Pind = 0.121; while for (i)-(l) the parameters are v/vf — 0.71, 
hy/ew = h z /e F = 0.1, fi/ep = 0.492, and P ind = 0.120. 
The blue dashed lines represent cuts of n s (k) along the 
k = (k x = 0,k y ,k z = 0) direction, while the red solid lines 
represents cuts of n s (k) along k = (k x , k y — 0, k z — 0). 



VII. DENSITY OF STATES 

The density of states for spin s can be written as 

p s (o;) = ^A s (k,o;), (28) 

k 

in terms of the spectral function A s (k, uj). An analy- 
sis of the spin-dependent density of states is useful to 
provide the frequency (energy) dependence of the spin- 
polarization of the system. 

In Fig. [51 we show the density of states p s (w) for vari- 
ous parameters with an energy broadening rj/ep = 0.01. 
We show specifically the case of zero spin-orbit coupling 
and Zeeman fields in Fig. [81(a), which has the charac- 
teristic square-root frequency dependence for a three- 
dimensional system. In this case the system is not polar- 
ized and the tails below the bottom of the energy band 
are due to the finite energy broadening. We show two 
other situations for comparison corresponding to cases 
which are polarized with finite Zeeman fields, as well as 
with a finite spin-orbit coupling. In Fig. [51(b)- (c), we 
show that the band edges shift to lower frequencies when 
the Zeeman and spin-orbit fields are turned on. Further- 
more, in Fig. [51(b) the spin-dependent density of states 
are shown for the case where parity is not violated, cor- 
responding to hy/ep = 0, h z /ep = 0.1, and v/vp = 0.71; 
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FIG. 7. (color online) Finite temperature (T/cf = 0.05) 
momentum distributions n-f(k) (two left-most columns) and 
ni(k) (two right-most columns). All panels corresponds to 
values of ERD spin-orbit coupling v/vf = 0.71, and Zeeman 
fields hy/eF = h z /eF = 0.1. For (a)-(d) the chemical poten- 
tial is /h/cf — 0.492, and the induced polarization is Pi n d = 
0.120. Similarly, for (e)-(h) p/e F = 0, and P ind = 0.143; while 
for (i)-(l) the parameters are /jl/cf = —0.37, and P in d = 0.105. 
The blue dashed lines represent cuts of n s (k) along the direc- 
tion k = (k x — 0,k y ,k z — 0) direction, while the red solid 
lines represents cuts of n s (k) along k = (k x , k y = 0, k z = 0). 



and in Fig. [5Jc) the spin-dependent density of states are 
shown for the case where parity is violated, corresponding 
to hy/ep = hy/ep = 0.1 and v/vp = 0.71. Since momen- 
tum is integrated over, there is no clear signature of par- 
ity violation in p s (uo) as there is in momentum-resolved 
observables such as the spectral density A s (\s.,lj), mo- 
mentum distribution n s (k), or helicity dispersions e^(k) 
and e^.(k) discussed earlier. However, we show the differ- 
ent spin-dependent density of states for comparison. The 
kinks present in these figures reflect the location of the 
maxima and minima of the helicity bands. 

Having discussed the density of states, we present next 
an analysis of thermodynamic properties. 
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FIG. 8. (color online) Spin-dependent density of states p s (uo) 
in units of N/cf for (a) v/vf = and h z /eF = h y /eF = 0, 
(b) v/vf = 0.71 and h y /cF = and h z /eF = 0.1, and (c) 
v/vf = 0.71 and h z /eF = h y /cF = 0.1. The blue dashed 
line corresponds to pt( CJ ) an d the solid red line corresponds 
to pi(uj). 



When no initial population imbalance is present, the total 
number of particles is fixed by 



N 



dp ; T y 



(31) 



which determines the chemical potential p. In addition, 
because an external Zeeman field is present, we can define 
the induced polarization 



ind 



where the number of particles in spin-state s is 
N s = / dconp(uj)A s (k, co). 



(32) 



(33) 



Next, we begin our analysis of thermodynamic vari- 
ables that could be measured using the techniques al- 
ready developed for ultra-cold fermions 0] in the ab- 
sence of spin-orbit coupling. In the discussion that fol- 
lows, we will cover the pressure, the chemical potential, 
the isothermal compressibility and the induced magneti- 
zation (spin-polarization). 



IX. PRESSURE AND ENTROPY 



VIII. THERMODYNAMIC PROPERTIES 

From the energy spectrum of the Hamiltonian Ho de- 
fined in Eq. (fT0]h we can obtain the partition function 

as 

Z = Tr [exp(-Ho/T)] , (29) 
and the corresponding thermodynamic potential 

tt = -TlnZ. (30) 



The pressure of ultra-cold fermions in the presence of 
spin-orbit coupling and crossed Zeeman fields is 

P = -^ = -|E ln ( 1 + e "^ (k)/T ). ( 34 ) 

k,Q! 

where ft is thermodynamic potential, and a = (ff , J|) is 
the helicity spin index. 

In Fig. [9j we show the scaled pressure PV/(Nep) as 
a function of the Zeeman fields h y /ep and h z /ep for 
fixed spin-orbit coupling v/vp- Notice that in the ab- 
sence of spin-orbit and Zeeman fields the pressure re- 
duces to the standard results of an non-interacting Fermi 
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gas PV/(Ncf) = 2/5. The pressure is an even func- 
tion of both h y and h z and is shown for a range of h y 
and h z varying from zero to 1.56^. The spin-orbit cou- 
pling for Figs. EJa)-(b) is v/vp = 0, and the pressure is 
completely isotropic in the h y -h z plane, since the helicity 
bands e a (k) are also isotropic in h y -h z plane and parity is 
preserved. However, in Fig. [9^c)-(d), where v/vf = 0.71, 
the pressure is anisotropic in the h y -h z plane, since the 
helicity bands e a (k) are now anisotropic in the h y -h z 
plane and are not invariant under parity. 
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FIG. 9. (color online) Pressure PV/Ncf as a function of 
hy/eF and h z /eF at finite temperature (T/cf = 0.05) for (a)- 
(b) v/vf = and (c)-(d) v/vf = 0.71. The red solid line 
represents PV/Ncf as a function of h y /cF at h z /eF = 0, 
and the blue dashed line represents PV/Ncf as a function of 
hz/eF at hy/eF — 0. 

We note in passing that the entropy S of the sys- 
tem can be easily extracted, by rewriting the thermody- 
namic potential as £1 = T ^ k In [1 — np(£ a (k))] , when 
expressed in terms of the Fermi function 7ip(£ a (k)). Us- 
ing the relation S = — (<9Q/<9T) y ^, leads to the final 
result (ks = 1) 

S =~ J2 W(m^(k))ln[n F (m^(k))]}, (35) 

ka,m=± 

which is nothing but the entropy of a non-interacting 
Fermi gas in the presence of spin-orbit coupling and 
crossed Zeeman fields. We will not show plots of the 
entropy, or of the specific heat, which can also be easily 
obtained, but rather discuss next the chemical potential 
and its dependence on the crossed Zeeman fields. 



particles given in Eq. (|3T|) . which can be rewritten as 
7V = ]Tn + (k), 



(36) 



where n+(k) is the total momentum distribution defined 
in Eq. (|26p . The behavior of jjj as a function of the Zee- 
man fields h y and h z is shown in Fig. [TUJ which uses the 
fact that fi is an even function of these variables. The 
range of the Zeeman fields is also from to 1.5^. The 
case where parity is preserved is shown if Fig. ITUTa) and 
(b), corresponding to v/vp = 0, such that the chemical 
potential is isotropic in the h y -h z plane. Similarly, in 
Fig. [TUTc) and (d), we show the case corresponding to 
v/vf = 0.71, where parity is violated for any finite Zee- 
man component h y . As a result the chemical potential 
ji/ep is anisotropic in h y -h z plane, because the helic- 
ity bands are neither even nor odd in momentum space. 
This demonstrates that the anisotropy of /i is a measure 
of parity violation. Another quantity that can be easily 
measured and that reveals a similar effect is the isother- 
mal compressibility to be discussed next. 
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FIG. 10. (color online) Finite temperature (T/ef 
chemical potential /jl/cf as a function of h z 
in (a)-(b) v/vf — and in (c)-(d) for v/vf = 0.71. The red 
solid line represents fi(hy,h z — 0)/cf and the blue dashed 
line represents ji(h y — 0,h z )/eF- 



XI. ISOTHERMAL COMPRESSIBILITY 



X. CHEMICAL POTENTIAL 

The chemical potential \i in the Grand-canonical en- 
semble is determined by fixing the average number of 



The isothermal compressibility can be obtained from 
the knowledge of the pressure, as is defined as kt = 
— (1 /V)(dP / 8V)t,n • But in the Grand-canonical ensem- 
ble, where we need to fix the average number of particles, 
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the isothermal compressibility can be directly written as 

—(f),,- 

Using the relation defined in Eq. ([36]) and notic- 
ing that the partial derivative {dnp{x)/ dx) T v = 
— (l/T)nF(x)riF(— x), is directly related to the Fermi 
function uf(x), it is possible to write the expression for 
the isothermal compressibility as 



X 2 *? 4E M&(k))n F (-&(k))] . 



(38) 



k,ai 



The experimental extraction of the isothermal com- 
pressibility from measurements of density fluctuations 
was suggested theoretically several years ago both in har- 
monically confined systems [35| and optical lattices [36| , 
and early improvements in the detection schemes of den- 
sity fluctuations [3TM39(] became sufficiently sensitive to 
extract this information from experimental data. In a re- 
cent experiment [40[ using laser speckles, the isothermal 
compressibility and the spin susceptibility were measured 
as a function of interaction parameter via the fluctuation 
dissipation theorem throughout the evolution from BCS 
to BEC superfluidity in balanced Fermi systems. 

The atomic compressibility can be measured using the 
fluctuation-dissipation theorem relating the fluctuation 
in the density of particles n = N/V, where TV = (N) 
is the average number of particles, and TV is particle- 
number operator. The relation between the isothermal 
compressibility kt and particle-number (density) fluctu- 
ations is given by the relation: 



V (TV 2 ) - (AQ 2 



(39) 



We see no major technical impediment to use tech- 
niques that are sensitive to spin-dependent density fluc- 
tuations in population imbalanced Fermi- Fermi mixtures 
with equal masses [4l|, |42| , whether the imbalance is cre- 
ated initially via radio-frequency fields or via artificial 
spin-orbit and Zeeman fields. Such analysis was shown 
to be theoretically possible even for Fermi mixtures with 
unequal masses |43|, |44|, and preliminary experimental 
results for these systems [SB], [46| seem to indicate that 
indeed the compressibility and spin susceptibility matrix 
elements can be directly extracted from the local density 
and density fluctuation profiles. 

Thus, we show in Fig. [11] the isothermal compressibil- 
ity kt as a function of both h y and h z for a range of h y 
and h z varying from zero to 1.56^, at fixed temperature 
(T = O.Obep)- The spin-orbit coupling for Figs. [TTl a)- 
(b) is v/vf = 0, and the compressibility is completely 
isotropic, as parity is preserved. However, in Fig. Qjjc)- 
(d), where vjvp = 0.71, parity is broken for h y ^ 0. 
This parity breaking is reflected in the helicity bands as 
discussed earlier, and manifests itself in the behavior of 



the compressibility versus (h y ,h z ) via the anisotropy re- 
vealed in Fig. HTT c)-(d). Again such anisotropy is a reflec- 
tion of the parity violation caused by the simultaneous 
presence of the ERD spin-orbit field and h y . Another 
important property that can be measured is the spin po- 
larization as a function of the Zeeman fields h yi h z for 
fixed spin-orbit coupling, which is discussed next. 
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FIG. 11. (color online) Finite temperature (T/ef — 0.05) 
isothermal compressibility N 2 kt in units of N/cf as a func- 
tion of h y and h z is shown in (a)-(b) for v/vf = and 
in (c)-(d) for v/vf = 0.71. The red solid line repre- 
sents N 2 KT(h y ,h z = 0) and the blue dashed line represents 
N 2 k,t (h y = 0,h z ). 



XII. INDUCED SPIN POLARIZATION 

We can also analyze the induced spin polarization 
in the presence of crossed Zeeman fields and spin-orbit 
coupling. The spin polarization along the i th direction 
(i = x, y, z) is given by the expectation value 



(40) 



Such general expression can be particularized for each 
component. For instance, the expectation value 

(s x ) = lJ2(( c U c ^ + ( c U c ^))- (4i) 



expressed in the original spin basis, can be rewritten in 
the helicity basis as 

(S x ) = 5>k(v k + v* k ) [n F (^(k)) - n F (&(k))] . 

(42) 
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Analogously the expectation value of the spin-operator 
S y in the original spin basis is 



(Sy) 



[< c kT c H> - ( c U c kt> 



(43) 



in the original spin basis can be written as 

(Sy) = ~\Y. U ^ ~ <) MCt(k)) - M&(k))] 

(44) 

in the helicity basis. Lastly, the expectation value of S z 
in the original spin basis is 

(S z ) = \j2((cUc kt )-(cl i c kl )), (45) 

k 

which can be rewritten in the helicity basis as 
(Sz) = \ J2« ~ H 2 ) K(^( k )) - • (46) 



Finally, the averages (S x ) and (S y ) can be expressed 
as real and imaginary parts of the transverse spin- 
polarization (S±) = (S x ) — i(S y ), defined this way to be 
compatible with the definition of ftj_(k) = h x {k)—ih y {k). 
The transverse spin-polarization takes the final form 



(Si.) 



2^|/i eff (k)| 



[n F (C t (k))-n F (^(k))]. (47) 



Correspondingly the longitudinal spin polarization can 
be written as 



(48) 



which is directly related to the induced population im- 
balance Pi n d by the expression 



ind 



2(&) 
N 



(49) 



where N = + N± is the total number of particles, as 
defined earlier. 

For ERD spin-orbit coupling with field Yierd = vk x y, 
the transverse field h±(k) has only the y-component. 
This means that (S x ) is identically zero for any value of 
h y and h z for any value of v ^ 0, given that h x {k) = 0. 
However, (S y ) is not identically zero for the ERD case 
above, unless parity is preserved in the helicity bands 
£ a (k), which means h y = 0. For any finite value of h y , 
(S y ) is non-zero. This behavior is revealed in Fig. [I2j 
where the expectation values (S y ) and (S z ) of spin po- 
larization are shown as a function of (h y /eF, h z /eF)- In 
particular, we show in Fig. [12] the finite temperature 
(T/ep = 0.05) induced spin polarizations per particle 
(S z )/N (two left-most columns) and (S y )/N (two right- 
most columns) for ERD spin-orbit parameter v/vp = 



from (a) through (d), and for v/vp = 0.71 from (e) 
through (h). In (a) and (e), the red solid line repre- 
sents (S z )/N as a function of h y jep at h z jep = 0, and 
the blue dashed line represents (S z )/N as a function of 
h z €F at h y /€F = 0. In (d) and (h), the red solid line 
represents (S y )/N as a function of hyjep at Ii z /cf = 0, 
and the blue dashed line represents (S y )/N as a function 
of h z /eF at hy/ep = 0. 

Having analyzed several thermodynamic properties for 
non-interacting quantum degenerate ultra-cold fermions, 
which already present some fundamental non-trivial 
properties such as the violation of parity, we discuss next 
the effects of interactions and how parity violation af- 
fects the pairing temperature of such fermions and the 
superfluid order parameter. 
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FIG. 12. (color online) Finite temperature (T/ef = 
0.05) induced spin polarizations per particle (S z )/N (left two 
columns) and (S y )/N (right two columns) are shown in (a)- 
(d) for v/vf = and (e)-(h) for v/vf — 0.71. In (a) and (e), 
the red solid line represents {S z )/N as a function of h y /cF at 
hz/eF = 0, and the blue dashed line represents (S z )/N as a 
function of h z /eF at h y /e = 0. In (d) and (h), the red solid 
line represents (S y )/N as a function of h y /eF at h z /eF = 0, 
and the blue dashed line represents (S y )/N as a function of 
h z /eF at h y /eF = 0. 



XIII. EFFECTS OF INTERACTIONS 

In this section, we analyze briefly the effects of inter- 
actions in the presence of spin-orbit and crossed Zeeman 
fields, with focus on the dependence of the pairing tem- 
perature with respect to the interaction parameter for 
given spin-orbit and Zeeman fields. We also present a 
short discussion about the effects of parity violation on 
the superfluid order parameter. 

The interaction Hamiltonian is 

Hi = -g J rfrV>j(r)Vl(r)V4(r)V t (r), (50) 

where g represents the strength of the contact interaction. 
Only s-wave scattering is considered in regards to the 
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original spin states t and |. Converting the interaction 
term into momentum space leads to 



% = -5E &t (i)H<i). 



(51) 



where the pair creation operator with center of mass 
momentum q is ^(q) = ^^{(k + ^/^V^ - k + q/2), 
and g can be expressed in terms of the scattering length 
through the Lippman-Schwinger relation 



V 



Vm 
47ra s 



(52) 



The interaction Hamiltonian Hi can be written in the 
helicity basis as 



ni = -g J2 B ^(q)^(q). 

qo;/37(5 



(53) 



where the indices a, /?, 7, £ cover ff and J| states. Pairing 
is now described by the operator 

£Mq) = ^A a/3 (k 1 ,k 2 )$ a (k 1 )$ /3 (k 2 ) (54) 



and its Hermitian conjugate, with momentum indices 
ki = k + q/2 and k2 = — k + q/2. The matrix 
A a ^(k + q/2, — k + q/2) is directly related to the ma- 
trix elements of the momentum dependent SU(2) rota- 
tion into the helicity basis, and reveals that the center of 
mass momentum ki+k2 = q and the relative momentum 
ki — k2 = 2k are coupled and no longer independent. 



XIV. TENSOR ORDER PARAMETER 

From Eq. (p3|) it is clear that pairing between fermions 
of momenta k x and k2 can occur within the same he- 
licity band (intra-helicity pairing) or between two differ- 
ent helicity bands (inter-helicity pairing). For pairing at 
q = 0, the order parameter for superfluidity is the tensor 
A a/3 (k) = A A^(k,-k), where A = -^E^B^O)), 
leading to components: 

A^(k) = A Ok^-k - ^k^-k) (55) 

for total helicity projection A = +1; 

A M (k) = -A (^k^-k + v*y*_ k ) 

(56) 

A M (k) = A (iik^-k + <^-k) 
for total helicity projection A = 0; and 

A M (k) = A (u^v\ - <^-k) (57) 



for total helicity projection A = — 1. 

It is very important to emphasize that for non-zero 
spin-orbit coupling and crossed Zeeman fields h y and h z , 



the order parameter tensor A a ^(k) does not have well 
defined parity. For instance, while A^(k) and A^(k) 
have odd parity, the matrix elements A^ (k) and A^ (k) 
do not have well defined parity. However, we may still de- 
fine singlet and triplet sectors for the helicity basis, such 
that the singlet sector As 5 o(k) = [A^(k) — A^(k)] /2 
has even parity and the triplet sector defined by the com- 
ponents A^(k), A T)0 (k) = [A^(k) + A^(k)] /2 and 
A^(k) have odd parity for any value of the ERD spin- 
orbit coupling v/vp and crossed Zeeman fields hy/cp and 
hz/ep. The preservation of parity in the singlet and 
triplet sectors is also true for the Rashba-only (RO) case, 
but the order parameter breaks time-reversal symmetry. 

Within the mean field approximation, the Hamiltonian 
matrix in the helicity basis is 



H MF (k) 











Ai, t (k) 



A tt (k) A M (k)\ 
A M (k) %(k) 

-Z-w J 



(58) 

This Hamiltonian matrix is traceless, therefore the sum of 
its eigenvalues is zero, however the eigenvalues of HMF(k) 
are not invariant under parity. By labeling the eigenval- 
ues as Ei(k), E^k), Es(k) and £^(k) in decreasing order 
of energy, and using the tracelessness condition then the 
sum E 3 (k) + E±(k) = - [#i(k) + ^(k)] but each eigen- 
value Ei(k) does not a well defined parity. Typically 
these eigenvalues are even in momentum space, but not 
here because the parity violation induced by the simulta- 
neous presence of the crossed Zeeman fields and the spin- 
orbit coupling, thus, in the present case Ei{— k) 7^ Ei(k). 
However a generalized particle-hole symmetry applies 
leading to E 2 (k) = -E 3 (-k) and £i(k) = -E 4 (-k). 

The eigenvalues in this parity violating case can be ob- 
tained analytically for any mixture of Rashba and Dres- 

selhaus terms from the determinant Det ujI — HMF(k) 
which leads to the characteristic quartic equation 

uj 4 + a 3 (k)cj 3 + a 2 {k)uj 2 + ai(k)cj + a (k) = (59) 

for each momentum k. Here, the coefficient of the cubic 
term 03 (k) = — Ei^(k) is the sum of the eigenvalues 
Ei(k) of the the Hamiltonian matrix HMF(k) and there- 
fore vanishes. The coefficient of the quadratic term is 
<22 (k) = J2i<j ^(k)£"j (k), while the coefficient of the lin- 
ear term is ai(k) = Yi^j^i ^(k^-^i^k). The last 
coefficient is just the product of the four eigenvalues lead- 
ing to a = £i(k)£ 2 (k)£ 3 (k)£ 4 (k). 

In the particular case of ERD spin-orbit coupling with 
crossed Zeeman fields, the coefficients become as = 0, 
and the coefficient of the quadratic term takes the form 



a 2 



2(^( k ) + |Ao| 2 + |^| 2 + |^| 2 + |^, 

(60) 

while the coefficient of the linear term is 



ai = -8K+(k)(vk x )h y , 



(61) 
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and lastly the coefficient of the zero-th order term is 
a = ^(k)^(k)^(-k)^(-k) + |A | 2 cg(k) (62) 

where ag(k) = {2K\ (k) + |A | 2 + ftg(k)) with ftg(k) = 

2|f ^| 2 — 2|/i 2y | 2 — 2|/i 2 | 2 . Here, if + has the same definition 
used in the paragraph that follows Eq. (|9j) , and is a mea- 
sure of the kinetic energy with respect to the chemical 
potential. 

Even in this simpler case of ERD spin-orbit coupling, 
the precise analytical form of the eigenvalues in the pres- 
ence of crossed Zeeman fields is quite cumbersome, and 
we do not list them here explicitly. Rather, we discuss 
next the consequences of parity violation on the pair- 
ing temperature of ultra-cold fermions in the presence of 
spin-orbit and crossed Zeeman fields. 



XV. PAIRING TEMPERATURE 

From the excitation spectrum discussed above, we ob- 
tain the corresponding thermodynamic potential as 



9 



k,j=l,4 



(63) 

from which the order parameter equation is determined 
via the minimization of Omf with respect to |Ao| 2 , lead- 
ing to 



V 1^ dEjjk) 



(64) 



where n F (Ej(k)) = l/(exp (Ej (k)/T) + 1) is the Fermi 
function for energy Ej(k). The contact interaction g can 
be elliminated in favor of the scattering length a s via the 
Lippman-Schwinger relation defined in Eq. ([52]) . 

The total number of particles N = + N± 

is defined from the thermodynamic relation TV = 
— (dftMF/dfi) T v , and leads to the corresponding num- 
ber equation 



(65) 



since the system is assumed to have no initial population 
imbalance. 

The self-consistent solutions of Eq. ([64]) and ([65]) guar- 
antee the existence of mean field solutions for the order 
parameter amplitude |Ao| and the chemical potential /i 
as a function of the Zeeman fields h y and h z , the spin- 
orbit coupling v and scattering length a s . However, the 
thermodynamic stability of the solutions obtained has 
to be tested against the maximum entropy condition (or 
minimum of the thermodynamic potential) over the same 
parameter space spanned by the variables h y , h z , v and 



a s , which determine the phase space of the present sys- 
tem. However, we discuss here only the effects of crossed 
Zeeman fields h y , h z on the pairing temperature T p ob- 
tained by solving the mean-field self-consistent relations 
defined by Eq. ([64]) and Eq. ([65]) with the order param- 
eter amplitude set to zero, i. e., Ao = 0. 

In Fig. [I3j we show the pairing temperature T p /ep as 
a function of ERD spin-orbit coupling v/vp for selected 
values of h y , h z and interaction parameter 1 / (kpa s ). The 
lines are guides to the eye given that the number of 
points does not form a dense set. In (a) h y /cp = 
and in (b) h y /ep = 0.08 both at l/(kpa s ) = —0.5, 
showing the behavior of T p on the BCS side of unitar- 
ity. In (c) hy/ep = and in (d) h y /ep = 0.4 both at 
l/(kpa s ) = 1-0, showing the behavior or T p on the BEC 
side of unitarity. In (a) and (b), the black-dotted line la- 
bels h z = 0, the red-dashed line h z jtp = 0.05, the green- 
dash-dotted line h z jep = 0.1, and the blue-solid line 
h z /ep = 0.15. However, in (c) and (d), the black-dotted 
line labels h z = 0, the red-dashed line h z /ep = 0.5, 
the green-dash-dotted line h z /ep = 1, and the blue-solid 
line h z /ep = 1.5. Notice that the relative suppression 
of pairing with zero center of mass momentum in the 
BCS side (l/(kpa s ) < 0) is larger than in the BEC side 
(l/(kpa s ) > 0), since it relies strongly on pairing of states 
only close to the Fermi wavevectors kp and — kp. 

In Fig. [HJ we show the pairing temperature T p /ep 
versus scattering parameter l/(kpa s ) for various values 
of the artificial Zeeman field components h y and h z and 
particular values of the ERD spin-orbit coupling v. The 
black cross line corresponds to parameters v/vp = 0, 
hy/cp = hz/ep = 0, where there are no ERD spin-orbit 
coupling and no Zeeman fields. The red solid line corre- 
sponds to parameters v/vp = 1, h y /ep = h z /ep = 0. No- 
tice that these lines coincide, because the ERD spin-orbit 
field can be gauged away producing exactly the same re- 
sults for any value oiv/vp so long as hy/cp = h z /ep = 0. 
The blue dashed-dotted line describes the case for pa- 
rameters v/vp = 1, hy/ep = 0, and h z /ep = 1, show- 
ing that the presence of Zeeman field h z in the BCS 
regime produces an energy cost for pairing of Fermions 
with opposite momenta and zero center-of-mass momen- 
tum, thus reducing the pairing temperature substantially. 
However, the purple dashed line describing the situation 
corresponding to v/vp = 1, h y /ep = 0.4, and h z /ep = 
shows a much stronger suppression of the pairing tem- 
perature than for the case of vjvp = 1, h y jep = 0, and 
h z /cF = 1 (blue dot-dashed line), because it becomes 
much more difficult for fermions with opposite momenta 
to pair with zero center-mass-momentum due to the par- 
ity violation in the excitation spectrum of the fermions 
introduced by h y when v/vp ^ 0. Lastly, the green 
dashed line shows the case of v/vp = 1, h y /ep = 0.4, 
and h z /ep = 1, where the combined effect of the Zee- 
man energy cost and parity violation lead to a dramatic 
reduction of the pairing temperature in the BCS region 
and even near unitarity. 

Another important point to emphasize in Fig. [TH re- 
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FIG. 13. (color online) Pairing temperature T p /cf as a 
function of ERD spin-orbit coupling v/vf for selected values 
of h y , h z and interaction parameter l/(kFCL s ). In (a) h y /cF = 
and in (b) h y /e F = 0.08 both for l/k F a s = -0.5. In (c) 
hy/eF = and in (d) h y /e F = 0.4 both for l/(k F a s ) = 1.0 For 
(a) and (b), the black dotted line labels h z — 0, the red dashed 
line h z /eF = 0.05, the green dash-dotted line h z /eF = 0.1, 
and the blue solid line h z /eF — 0.15. In contrast, for (c) and 
(d), the black dotted line labels h z = 0, the red dashed line 
h z /eF = 0.5, the green dash-dotted line h z /eF = 1, and the 
blue solid line h z /eF — 1.5. Notice that it is much easier to 
suppress pairing in the BCS side (l/(kFa s ) < 0), than in the 
BEC side (l/(k F a s ) > 0). 



lates to the bending of T p in the two last cases corre- 
sponding to the dashed purple and dashed green lines. 
This bending indicate and instability of zero center-of- 
mass momentum pairing towards finite center-of-mass 
momentum pairing, and point of infinite slope in both 
curves indicates the separation of the two regimes. In 
the cases where the pairing temperature reflects the crit- 
ical temperature of the system, the locations of infinite 
slope would correspond to a Lifshitz point, and a small 
region to the left of the negative slope regime would cor- 
respond to a superfluid with finite center of mass mo- 
mentum, which is favored due to parity violation in the 
helicity bands. 

For the parameters discussed the pairing temperature 
is not largely affected in the BEC regime, since the bind- 
ing of fermions is controlled by the emergence of two- 
body bound states with binding energy E bl and not by 
Cooper pairing in the presence of a Fermi sea. In the 
BEC regime, an estimate of T p can be given by con- 
sidering the chemical equilibrium condition /i# = 
where \±b is the chemical potential of the bosons formed 
by tightly bound fermions, and \i F is the chemical po- 
tential of unbound fermions. Since T p /e F > 1, both 



bosons and unbound fermions are highly non-degenerate, 
and behave like classical ideal gases, in which case the 
chemical potential can be directly calculated and used 
to obtain the relation T p /e F ~ \E b /e F \/ \n\E b /e F \ 3 ^ 2 to 
logarithmic accuracy. Here, the binding energy E b /e F 
is a function of the interaction parameter l/(k F a s ), the 
ERD spin-orbit parameter v/v F and the Zeeman fields 
h y /e F and h z /e F . The logarithmic term is an en- 
tropy correction that reduces the pairing temperature 
T p to a value much lower than the absolute value of 
the binding energy E b of two fermions. The pairing 
temperature T p /e F is essentially the same for values of 
l/(k F a s ) > 2 (not shown in Fig. [12]), because the bind- 
ing energy \E b /e F \ ^> max.[v/vF, h y /e Fl h z /e F ], in which 
case \E b /e F \ « 2/{k F a s ) 2 and the pairing temperature 
is T p /e F ~ [2/ (k F a s ) 2 )/ ln[l/(k F a s )] 3 , which agrees with 
the numerical calculations in the regime of l/(k F a s ) ^> 1. 

We will not present here a discussion of the conden- 
sation or critical temperature T c as it requires a full 
calculation of fluctuations effects, which is now under- 
way [47l]. However, in the extreme BEC regime, where 
l/(k F a s ) ^> 1, the critical temperature T c can be ob- 
tained from the Bose-Einstein condensation temperature 
T c « Tbec = [n B /C(3/2)] 2 / 3 27r/ra B , where the zeta 
function C(3/2) « 2.6124, the density of bosons tlb = n/2 
is half the density n of fermions, and is the mass of 
bosons which is a function of the mass m of fermions, and 
the parameters u, h y and h z . The ratio m^/m can be pa- 
rameterized by the dimensionless ratios v/v Fl h y /e F and 
h z /e F . In the limit of l/(k F a s ) — ^ 00 with finite v/v F , 
h y /e F and h z /e Fl the Boson mass 4 2m, but cor- 
rections depending on the aforementioned ratios tend to 
make the mass heavier, and for finite but large l/(k F a s ) 
this mass increase tends to reduce Tbec- However, the 
effect of interactions between the effective bosons is very 
subtle and the understanding of its dependence on the 
dimensionless ratios l/(k F a s ), v/y F , h y /e Fl and h z /e F 
is being currently investigated [47j . 

After a detailed discussion of the effects of parity vio- 
lation on interactions, order parameter, and pairing tem- 
perature T p , which is is dramatically affected in the BCS 
regime, but largely unaffected in the BEC regime for the 
parameter range investigated, we are ready to present 
our conclusions. 



XVI. CONCLUSIONS 

We have analyzed the normal state of a degenerate 
Fermi gas in the presence of artificial spin-orbit cou- 
pling and crossed Zeeman fields. The specific form of the 
spin-orbit field chosen corresponds to a mixture of equal 
Rashba and Dresselhaus (ERD) terms h^^(k) = vk x y, 
which has been experimentally realized. The artificial 
Zeeman field h z along a defined spin quantization z cor- 
responds to the Raman intensity of the laser beams, and 
the crossed Zeeman field h y pointing along the same di- 
rection the ERD spin-orbit field corresponds to the fre- 
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1/kF a s 

FIG. 14. (color online) The plot of the pairing tempera- 
ture Tp/eF as a function of l/kF(is> The black dotted line 
represent the case for ERD spin-orbit coupling v/vf = 0, 
and Zeeman fields h y /cF — h z /cF = 0; the red solid line 
for corresponds to v/vf = 1, h y /cF = 0, and h z /eF = 0; 
the blue dashed-dotted line to v/vf = 1, h y /cF = 0, and 
hz/eF = 1; the purple dashed line describes the case of 
v/vf = 1, hy/eF = 0.4, and h z = 0; and finally the green 
dashed line corresponds to v/vf = 1, h y /eF = 0.4, and 
h z /eF - 1. 

quency detuning from the atomic transition coupling two 
spin states. 

In such configuration, the eigenvalues of the non- 
interacting problem are obtained in a generalized helic- 
ity basis, where the helicity spin projection points ei- 
ther along or opposite to the quantization axis denned 
by the effective magnetic field h e ff(k) = (0,h y + vk x , h z ). 
In this case, we have shown that the presence of spin- 
orbit and crossed Zeeman fields lead to parity violation 
in the excitation spectrum of the non-interacting Fermi 
gas, which has immediate consequences for the Fermi sur- 



face, spectral density and momentum distribution, which 
also do not have well defined parity. Such parity viola- 
tion emerges in momentum-resolved spectroscopic quan- 
tities because they are all functions of the magnitude 
|h e ff(k)| of the effective field h e ff(k) = (0,h y + vk x ,h z ), 
which does not have well defined parity, since h e ff (— k) = 
(0,h y - vk x ,h z ). 

In addition, we have shown that information on par- 
ity violation can be extracted from momentum- averaged 
thermodynamic properties such as pressure, entropy, 
chemical potential, compressibility, and spin-polarization 
as a function of crossed Zeeman field components hy/ep 
and h z /ep for fixed spin-orbit coupling v/vp &t fixed tem- 
perature T. A signature of parity violation is the differ- 
ent behavior of any chosen thermodynamic quantity as 
a function of h y for h z = and as a function of h z for 
h y = 0. This anisotropy of thermodynamic properties in 
the (h y , h z ) plane for finite spin-orbit coupling v/vp is a 
direct reflection of the lack of parity (inversion symme- 
try). 

Lastly, we analyzed the effects of interactions for a de- 
generate Fermi gas when parity is broken, and investi- 
gated how parity violation influences the order parame- 
ter of a uniform superfluid and the fermion pairing tem- 
perature. We noticed that the order parameter tensor 
for a uniform superfluid in the generalized helicity basis 
no longer possesses inversion symmetry, however singlet 
and triplet pairing in can still be defined this basis and 
preserves parity. Furthermore, we found that the effects 
of parity violation are strong in the pairing temperature 
of fermions, because it becomes increasingly more diffi- 
cult to pair states with zero center of mass momentum 
between helicity bands with progressively larger loss of 
inversion symmetry. 
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